Prediction of stable insulating intermetallic compounds 
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We explore the stability of structure exhibiting hybridization gaps across a broad range of binary 
and ternary intermetallic compositions by means of band structure and total energy calculations. 
This search reveals previously unknown metal-based insulators, some with large gaps exceeding 1 
eV, such as Ab>Fe and AUIrRe. We confirm large gaps using a hybrid density functional including 
exact exchange, and predict a gap of 2.2 eV for AlMnSi in the Pearson type tP6 structure, which is 
a chemically ordered ternary variant of the prototype MoSi2 (Pearson type tI6) structure. 

PACS numbers: 



I. INTRODUCTION 

Certain intermetallic alloys are electrically insulating 
despite the good metallic characters of their constituent 
elements. A classic example of this is provided by the 
compound AI2RU which is insulating as result of gaps in 
its density of states around its Fermi level Such ma- 
terials pose intriguing questions as to the origins of their 
insulating properties and additionally may find potential 
applications, for example as thermoelectric materials |2|. 

Previous studies 0-0| have explained the origins of 
this gap, and generalized the families of compounds in 
which it is expected to occur. The essential mechanism 
is to establish a gap between bonding and antibonding 
states through the creation and occupation of hybridized 
orbitals, then to exactly fill the bonding states. For struc- 
tures based on the prototype MoSi2 (Pearson type tI6) 
and TiSi2 (Pearson type oF24) this occurs at composition 
(Al) 2 (Fe), where the parenthesis around (Al) indicates an 
aluminum-group metal of valence three (specifically Al or 
Ga) and (Fe) is an iron-group transition metal of valence 
eight (i.e. Fe, Ru or Os). Notice these choices result in 
a total of fourteen valence electrons per transition metal 
atom. 

Gap width may be controlled by selection of the spe- 
cific combinations of elements, or by alloying with adja- 
cent elements in the periodic table (e.g. replacing the 
(Fe)-group element with a combination of (Mn)-group 
and (Co)-group), while maintaining constant valence per 
transition metal atom. Alternatively, we may replace half 
of the (Al)-group elements with (Si)-group, while simul- 
taneously replacing the (Fe)-group transition metal with 
a (Mn)-group element. 

Here we report a thorough study of the tI6 and oF24 
alloy families. The crucial new feature of our work is 
we simultaneously evaluate the compounds' enthalpies of 
formation, so that we may predict cases where the pro- 
posed structures are likely to form as stable compounds, 
as bandgap engineering through chemical substitution re- 
quires monitoring both variation of band gaps as well as 



thermodynamic stability relative to competing phase for- 
mation [8]. In this manner, we predict the occurrence of 
thirteen previously unknown compounds, each of which 
exhibit electronic band gaps. We then checked the gaps 
predicted by conventional DFT utilizing a hybrid func- 
tional including exact exchange. Among our newly pre- 
dicted stable compounds is Al-Fe, with a gap of 1.1 eV, 
Al-Ir-Re, with a gap of 1.3 eV, and Al-Mn-Si with a gap 
of 2.2 eV. 



II. STRUCTURES 

AI2RU has been reported to occur in both the tI6 Q 
and the oF24 [n| structures. Indeed, both structures 
share common structural elements, and are stabilized by 
similar mechanisms. The essential structural feature is 
an Ru-centered ring of six Al atoms (see Fig. [T|), re- 
peated periodically to tile the plane. A third structure, 
based on the CrSi2 prototype (hP9) shares a similar mo- 
tif. The ring is a slightly irregular hexagon, with only 
a 2-fold rotational symmetry. This ring lies perpendic- 
ular to the [110] axis in tI6 and the [001] axis in hP9 
and oF24. The differences between structures rests in 
the stacking of these TM-centered rings (TM stands for 
transition metal). Consider a triangle of three TM atoms 
within a layer. There are four possible sites for the TM 
in an adjacent layer, either the vertex site or else one of 
the three edge centers of this triangle. In tI6 the TM 
alternates between two such sites, in hP9 between three 
and in oF24 all four sites are utilized. 

Chemical ordering in ternary compounds lowers the 
symmetry from its binary prototype. For example, 
the MoSi2 prototype is Pearson type tI6, space group 
I4/mmm, while the optimal decoration for AlMnSi is 
Pearson type tP6 (Primitive, not body centered) with 
space group P4/nmm . The chemical occupation is uni- 
form within layers when the structure is viewed down 
the 4-fold axis. For example, an AJ^Mn layer alternates 
with an Si2Mn layer (see Fig. [IJ)), rather than mixing Al 
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FIG. 1: Structures of the class. Colors indicate atomic 
species, size indicates vertical height, (a) AI2OS in the Pear- 
son tI6 (MoSi2) structure. TM hexagons in alternate layers 
outlined in solid and dashed lines, (b) AlMnSi in tP6 struc- 
ture. View along (001) axis reveals uniform chemical identity 
of layers, and relationship to cP2 (CsCl) structure, (c) CrSi2 
in the hP9 structure revealing three stacked hexagons, (d) 
Generic hexagon motif in single layer, (e) AI2RU in oF24 
(Si2Ti) structure reveals four stacked hexagon motifs, (f) 
Ternary decoration of AUIrRe in oF24 structure. 



and Si within a layer. Likewise the Si2Ti prototype of 
Pearson type oF24 has space group Fddd while the opti- 
mal decoration for AlIrRe maintains Pearson type oF24 
but lowers the space group to F222 (Fig[T]/). Energeti- 
cally optimized Wyckoff positions for the two prototypi- 
cal structures are given in Table HI 



III. METHODS 

We utilize VASP [ll]> EH to carry out first principles 
total energy calculations. All atomic positions and lat- 
tice parameters are fully relaxed. We increase our k-point 
densities until energies have converged to 1 meV/atom. 
Default energy cutoffs are employed, which are sufficient 
to converge energy differences to this tolerance. Because 
of the volume relaxation, our energy may be considered 



as an enthalpy at P=0 and T=0K. We adopt projector 
augmented wave potentials [H, H3 |. and for a density 
functional we choose the PBE generalized gradient ap- 
proximation. Spin polarization is utilized in appropriate 
cases. 

Conventional density functional theory generally un- 
derestimates band gap s, so we employ a hybrid func- 
tional [15] . HSE06 [16}, that includes a portion of ex- 
act exchange, to obtain more reliable values. However, 
the hybrid functional is computationally expensive so we 
do this only in selected cases. We also report pseu- 
dogap widths, which we define as the energy interval 
within which the density of states remains below 0.20 
states/ev/atom. 

Our stability calculations follow methods outlined 
in (l7j . Primarily, we calculate enthalpy of formation 
AH by subtracting the calculated enthalpies from the 
enthalpies of the pure elements in their stable forms. In 
addition to the enthalpies of our tI6 and oF24 structures, 
we include all known structures within the given alloy 
system as well as selected hypothetical structures that 
occur in chemically similar systems. 

In the case of ternary systems, the atomic decoration is 
not constrained, due to the mixing of TM species, or sub- 
stituting (Si)-group for (Al)-group elements. To search 
for optimal patterns of chemical ordering, we first ex- 
haustively enumerate all symmetry-independent configu- 
rations within cells specified below. We examine for two 
representative alloy systems of each class, to check for 
consistency of the optimal chemical decoration, and then 
subsequently apply the optimal decoration to all other 
ternary alloy systems from the corresponding family. 

(B/C)Mn class. Our representative alloy systems were 
AlMnSi and AlReSi. In the MoSi 2 .tI6 prototype struc- 
ture, we limited our ground state search to 1x1x2 and 
\[2 x \/2x 1 supercells. Each of these had 10 independent 
configurations. In the Si2Ti.oF24 prototype case, the 
unit cell with 8 Si and 8 Al atoms occupying "Si" sites 
yields 293 independent configurations. We limited our 
ground-state search to configurations of relatively high 
symmetry. 

Al(Mn/Co) class. Our representative alloy systems 
were AlIrRe and AlReRh. In the MoSi2.tI6 prototype 
structure, there are 2 independent configurations each in 
the 1x1x2 and %/2x ^2x1 supercells. In the Si 2 Ti.oF24 
prototype case, there were 8 independent assignments for 
transition metals on the "Ti" sites. 

Given a database of relaxed structures, we calculate 
the convex hull of enthalpy as a function of composi- 
tion. Predicted stable structures occupy the vertices of 
the convex hull. For structures that are unstable, we re- 
port the energy AE by which they lie above the convex 
hull. For stable compounds, negative values of AE indi- 
cate the formation enthalpy with respect to compounds 
of adjacent compositions. 

Fig. [5] gives an example for the case of Al-Mn-Si. We 
reproduce the stability 18] of most known low temper- 
ature binaries, and two of the ternaries. One of these 
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FIG. 2: (color online) Phase diagram of Al-Mn-Si. Heavy cir- 
cles indicate known stable phases, light circles indicate known 
high temperature phases, squares indicate unknown struc- 
tures or compositions. AH values are enthalpy of formation 
(meV/atom) of stable structures. Positive AE values are en- 
ergies above the convex hull for unstable structures, while 
negative AE is the formation enthalpy with respect to adja- 
cent compositions for stable structures. 
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Al 


2c 0.75 


0.75 


0.9014 


Ir 


4a 0.0 


0.0 0.0 


Mn 


2c 0.75 


0.75 


0.2736 


Re 


4d 0.75 


0.75 0.75 


Si 


2c 0.75 


0.75 


0.5846 


All 


8f 0.50 


0.6706 0.50 










A12 


8i 0.25 


0.0817 0.25 



TABLE I: Atomic structures of AlMnSi.tP6 (left block; 
space group P4/nmm (#129), a=3.027A, c=7.932A) and 
Al 4 ReIr.oF24 (right block; group F222 (#22), a=4.745A, 
6=8.139A, c=8.925A). Wyckoff sites are labeled by site mul- 
tiplicities. 



IV. RESULTS 

Table [TTI presents our predicted stable structures along 
with their band gaps (if any) and pseudogap widths. 
Complete data is given in Table [TTJ We confirm the sta- 
bility of A12Ru, A120s, Ga2Ru and Ga20s in the oF24 
structure, and AlGeRe in the tI6 structure. We pre- 
dict the stability of (Al,Ga)4TcRh, (Al,Ga)4ReIr and 
AlGeMn in the oF24 structure and A12Fe, AlSiMn, 
(Al,Ga)SiTc, GaSiRe, (Al,Ga)GeTc and GaGeRe in the 
tI6 structure. 

We previously discussed the case of A12Fe [19]. Al- 
though total energy calculations predict it to be stable, 
the observed A^Fe phase is intrinsically disordered with 
low symmetry. This phase may be thermodynamically 
stabilized at high temperature by its vibrational entropy. 
It is probable that formation of the true low temperature 
stable phase is kinetically hindered at low temperatures. 
The most interesting new prediction may be AlMnSi, as 
it is comprised of common elements and it has the largest 
predicted true gap (2.2 eV, 1 eV larger than elemental 
silicon) 

The electronic density of states and band structure of 
our predicted AlMnSi structure are shown in Fig. [3] As 
expected, the conventional DFT gap is severely under- 
estimated, owing to the derivative discontinuity and de- 
localization error [2(|. Inspecting the band structure, 
we see that the band gap is indirect. Projecting onto 
atomic orbitals we find the top of the valence band at 
r is purely d xy localized on the Mn atoms, while the 
bottom of the conduction band, located along r — X is a 
hybrid of Al and Si p x and p y with Mn d xz and d yz states. 
A nearly degenerate conduction band local minimum of 
similar hybridized character occurs along M — T. An- 
other interesting feature is the nearly fiat valence band 
minimum along T — Z , which is responsible for the step 
function-like onset of the density of states near E = —15. 



V. DISCUSSION 



ternaries (Ali6Mn4Si3.cP138) is an icosahedral quasicrys- 
tal approximant, and exhibits a gap in its density of 
states. However, two reported phases (Alo.6Mno.sSio.4 
and AlgM^Sig) have unknown structure and hence are 
omitted from our study, and three reported phases (Pear- 
son types cP8, oF24 and hP9) exhibit partial site occu- 
pancy or mixed chemical occupancy and lie off their ideal 
stoichiometries, and hence are presumed to be stable only 
at high temperatures. Meanwhile, our predicted lowest 
energy of AlMnSi in the tP6 structure has not been re- 
ported experimentally. We suspect the prediction is cor- 
rect, as the AE — — 48 meV/atom is well beyond the 
expected uncertainty of DFT, and our exact exchange 
calculation confirms stability of tP6. Perhaps oF24 is a 
high temperature phase whose entropy inhibits formation 
of the true low temperature tP6 phase through conven- 
tional sample preparation methods. 



We have predicted a large number of stable intermetal- 
lic compounds with hybridization gaps based on chemical 
ordering on underlying tI6 and oF24 structures, several 
of which have not been previously reported. In some 
cases they may not have been previously reported be- 
cause the alloy systems have not been throughly evalu- 
ated. In the case of Al-Mn-Si, a chemically disordered 
phase of structure type oF24 is observed, while we iden- 
tify the stable state as tP6 (a chemically ordered ternary 
variant of tI6). Presumably the entropy of chemical sub- 
stitution and atomic vibrations stabilizes the disordered 
structure at high temperature, making identification of 
the low temperature stable structure difficult. Prelimi- 
nary investigation confirms that both configurational and 
vibrational entropy of oF24 exceed that of tP6. While 
the oF24 structure also exhibits a gap, it is smaller than 
that of tP6, and the chemical disorder can be expected to 
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FIG. 3: (a) Electronic density of states of AlMnSi in the 
MoSi 2 .tI6 structure. Shown are DOS for the PBE and HSE06 
functionals with Ef = 0. Threshold of 0.2 states/eV/atom 
defines the pseudogap. (b) Band structure of AlMnSi in the 
MoSi2.tI6 structure calculated using the HSE06 functional. 
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Ab(Fe) 


PBE HSE06 


Ga 2 (Fe) 


PBE HSE06 


Chem Struct 


A_£/ A// Eg Ep Eg E p 


Chem Struct 


A.£/ A_£^ Eg Ep Eg Ep 


Al-Fc ° rth ° 
A1 * e Tetra 


22 -335 0.220 0.650 0.739 1.287 
-21 -358 0.008 0.401 1.115 1.663 


Ga-Fe ° rth ° 
^ a * e Tetra 


127 -94 0.141 0.469 0.614 1.115 
182 -39 0.011 0.016 0.001 0.016 


Al-Ru ° Tth0 
Tetra 


-17 -685 0.112 1.119 0.591 1.669 
33 -653 0.005 0.556 0.489 1.319 


r> Ortho 

Ga-Ru _ 

Tetra 


-20 -377 0.107 0.928 0.559 1.448 
74 -303 0.014 0.527 0.274 1.149 


Al-Os ° r t h ° 
Tetra 


-0.7 -560 0.698 1.287 1.061 2.060 
0.7 -559 0.354 1.151 1.057 1.882 


„ „ Ortho 
Ga-Os _ 

Tetra 


76 -156 0.644 1.142 1.199 1.767 
107 -125 0.008 0.987 0.669 1.752 


Al 4 (Mn)(Co) 


PBE HSE06 


Ga 4 (Mn)(Co) 


PBE HSE06 


Chem Struct 


Ai? A// Eg Ep Eg Ep 


Chem Struct 


A£/ A_ff Eg Ep Eg Ep 


Al-CoMn ° rth ° 
Al OoMn Tetra 


43 -364 0.274 0.700 0.690 1.413 
39 -368 0.028 0.394 0.605 1.545 


Ga-CoMn 0rth ° 
Ietra 


93 -144 0.277 0.554 x x 
163 -74 0.014 0.241 x x 


Al-TcRh 2 r * h0 
Tetra 


-50 -693 0.425 1.125 x x 
50 -643 0.013 0.528 0.023 1.199 


Ga-TcRh 9? ho 
Tetra 


-76 -402 0.307 1.024 0.805 1.666 
76 -326 0.014 0.429 x x 


Al-Relr 2 r f th0 
Tetra 


-40 -597 0.702 1.326 1.260 1.999 
85 -513 0.000 0.009 0.001 0.522 


Ortho 

Ga-Relr _, , 

Tetra 


-18 -234 0.905 1.250 1.426 1.833 
57 -177 0.009 0.751 x x 


Al(Si)(Mn) 


PBE HSE06 


Ga(Si)(Mn) 


PBE HSE06 


Aim Mr, ° rth ° 
Tetra 


50 -356 0.48 0.84 1.195 1.353 
-48 -403.0 0.435 0.821 2.156 2.254 


P„oi M „ Ortho 
Ietra 


134 -169 0.599 0.614 x x 
54 -249 0.370 0.681 1.761 2.022 


Aiq: T _ 0rth0 

AlSi-Tc Tetra 


76 -528 1.101 1.174 1.757 1.829 
-76 -604 0.485 1.351 1.320 2.498 


p.c: T . 0rtno 

LraOl-lC . 

Tetra 


30 -382 0.925 1.025 x x 
-30 -412 0.363 1.194 1.181 2.223 


Aim Re ° rth ° 
Ietra 


72 -395 0.336 0.954 1.050 1.614 
-33 -467 0.277 1.615 0.862 2.637 


n.Ql d„ 0rtno 
GaSi-Re Tetra 


132 -164 0.323 0.856 0.902 1.483 
-54 -296 0.135 1.230 0.908 2.295 


AlGe-Mn ° rtho 
Tetra 


-20 -164 0.612 0.645 0.000 0.000 
20 -144 0.478 0.638 x x 


GaGe-Mn ° rtho 
Tetra 


90 26 0.504 0.541 x x 
68 3 0.212 0.489 x x 


AlGe-Tc J 1 ? 
Tetra 


21 -402 0.754 0.953 x x 
-21 -423 0.311 1.089 1.151 2.107 


GaGe-Tc 2 rt + ho 
Tetra 


19 -216 0.694 0.872 x x 
-19 -235 0.178 0.869 0.880 1.816 


r> Ortho 
AlGe-Re ^ 


63 -178 0.172 0.781 x x 
-43 -241 0.106 1.415 0.870 2.284 


n n -d Ortho 
GaGe-Re ^ , 

Tetra 


71 40 0.199 0.745 x x 
-20 -30 0.019 0.861 0.811 1.907 



TABLE II: Stability and bandgap data for various alloy systems. Element names in parenthesis indicate columns of the periodic 
table. AH and AE are enthalpy of formation and instability energy in units of meV/atom. E g and E p are the energy gap 
and pseudogap in units of eV. PBE and HSE06 are the density functionals. Newly predicted stable structures are in bold. 
Previously known and reconfirmed stable structures are in italics. Maximum gap and pseudogap are in bold. 



